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Abstract. In this note we survey some recent results for the Euler 
equations in compressible and incompressible fluid dynamics. The main 
point of all these theorems is the surprising fact that a suitable variant 
of Gromov's /i-principle holds in several cases. 



1. Introduction 

The starting point of this note is the Cauchy problem for the incompress- 
ible Euler equations, 

dtv + d\v x v ® v + Vp = 0, 

div^ v = 0, (1) 
v(0, •) = v , 

where the unknowns v and p are, respectively, a vectorfield and a scalar 
function defined on M n x [0, T). These fundamental equations were derived 
over 250 years ago by Euler and since then have played a major role in 
fluid dynamics. There are several oustanding open problems connected to 
([1]) and an extensive literature about them (see for instance the surveys [3], 
|23| , [39] ) . In three space dimensions little is known about smooth solutions 
apart from classical short time existence and uniqueness (see \50\ p 422]) and 
the celebrated blow-up criterion of Beale-Kato-Majda [4] (further references 
for both results are [35] , [2] , [17J and also the books [52] and [53] ) . On the 
other hand weak solutions are known to be badly behaved from the point 
of view of Hadamard's well-posedness theory: in the groundbreaking paper 
|66j V. Scheffer proved the existence of a nontrivial weak solution compactly 
supported in time. Nevertheless, weak solutions have been studied for their 
expected relevance to the theory of turbulence (see [60], [21], [37], [69]). 

In this survey we argue that the nonuniqueness of weak solutions for 
the incompressible Euler equations should be viewed as a suitable variant 
of the /i-principle. The original /i-principle of Gromov pertains to various 
problems in differential geometry, where one expects high flexibility of the 
moduli space of solutions due to the underdetermined nature of the problem. 
It was not expected that the same principle and similar methods could be 
applied to problems in mathematical physics (we quote Gromov's speech 
at the Balzan Prize [41]: The class of infinitesimal laws subjugated by the 
homotopy principle is wide, but it does not include most partial differential 
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equations ( expressing infinitesimal laws ) of physics with a few exceptions in 
favour of this principle leading to unexpected solutions. In fact, the presence 
of the h-principle would invalidate the very idea of a physical law as it yields 
very limited global information effected by the infinitesimal data. See also 
the introduction in the book |36j). 

As pointed out in the important paper [57] by S. Miiller and V. Sverak 
the existence theory for differential inclusions (see also \17\ fT8j [271 H5| [4"6] ) 
has a close relation to the ^-principle. In particular the method of convex 
integration, introduced by Gromov and extended by Miiller and Sverak to 
Lipschitz mappings, provides a very powerful tool to construct solutions. In 
the paper [30j these tools were suitably modifed and used for the first time to 
explain Scheffer's nonuniqueness theorem. It was also noticed immediately 
that this approach allows to go way beyond the result of Scheffer. Indeed it 
has lead to new developments for several equations in fluid dynamics. The 
goals of this note is to survey these results, list some open questions and 
point our some new directions. The most exciting ones address the question 
of whether these methods might give us some better understanding of the 
theory of turbulence. 

2. Weak solutions 

We start with a survey of the non-uniqueness theorems for weak solutions 
of (pQ). By a weak solution we mean, as usual, an Lf oc vectorfield which 
solves the equations in the sense of distributions. 

Definition 2.1. A vectorfield v G L 2 oc (R n x (0,T)) is a weak solution of the 
incompressible Euler equations if 

I [ d t if ■ v + Vip : v <g> v dxdt = (2) 

Jo JR n 

for all p G C c °°(lfr x (0, T); R n ) with div (p = and 

r-T 

v ■ Vip dxdt = for all ip G C 6 °° (R n x (0, T)). (3) 



o 



When vq G L 2 (M n ), the vectorfield v is a weak solution of ([T]) if ([2]) can be 
replaced by 



/ / dtcp ■ v + V</2 : v ® v dxdt + / ip(x, 0) • vq{x) dx = (4) 
Jo ii™ JK. n 

for all ip G C c °°(M n x [0,T);]R n ) with div ip = 0. 

2.1. Weak solutions with compact support in time. As already men- 
tioned, the first nonuniqueness result for weak solutions of (pQ) is due to 
Scheffer in his groundbreaking paper [66] . The main theorem of |66| states 
the existence of a nontrivial weak solution in L 2 (M 2 x R) with compact sup- 
port in space and time. Later on Shnirelman in [68] gave a different proof of 
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the existence of a nontrivial weak solution in L 2 (T 2 x M) with compact sup- 
port in time. In these constructions it is not clear if the solution belongs to 
the energy space (we refer to the next section for the relevant definition). In 
the paper [30J wc provided a relatively simple proof of the following stronger 
statement. 

Theorem 2.2 (Non-uniqueness of weak solutions). There exist infinitely 
many compactly supported weak solutions of the incompressible Euler equa- 
tions in any space dimension. In particular there are infinitely many solu- 
tions v G L°° n I? fo (p} for v o = and arbitrary n > 2. 

2.2. The Reynolds stress and subsolutions. The proof in [30] is based 
on the notion of subsolution. In order to motivate its definition, let us recall 
the concept of Reynolds stress. It is generally accepted that the appearance 
of high-frequency oscillations in the velocity field is the main responsible for 
turbulent phenomena in incompressible flows. One related major problem is 
therefore to understand the dynamics of the coarse-grained, in other words 
macroscopically averaged, velocity field. If v denotes the macroscopically 
averaged velocity field, then it satisfies 

d t v + div (v ® v + R) + Vp = , , 

divT7=0, [) 

where 

R = v ® v — v®v. 

The latter quantity is called Reynolds stress and arises because the aver- 
aging does not commute with the nonlinearity v ® v. On this formal level 
the precise definition of averaging plays no role, be it long-time averages, 
ensemble- averages or local space-time averages. The latter can be inter- 
preted as taking weak limits and indeed,. Indeed weak limits of Leray so- 
lutions of the Navier-Stokes equations with vanishing viscosity have been 
proposed in the literature as a deterministic approach to turbulence (see 

B21, EH 0, El)- 

A slightly more general version of this type of averaging follows the frame- 
work introduced by Tartar [781 179] and DiPerna [35] in the context of conser- 
vation laws. We start by separating the linear equations from the nonlinear 
constitutive relations. Accordingly, we write ([5]) as 

dt v + div u + Vq = 
divv = 0, 

where u is the traceless part oiv®v + R. Since one can write 

R = (v — v) ® (v — v), 

it is clear that R > 0, i.e. R is a symmetric positive semidefinite matrix. In 
terms of the coarse-grained variables (v, u) this inequality can be written as 

v ® v — u < §e J, 
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where / is the n x n identity matrix and 



— 1 I 12 

is the macroscopic kinetic energy density. Motivated by these calculations, 
we define subsolutions as follows. Since they will appear often, we introduce 
the notation <Sq xn for the vector space of symmetric traceless nxn matrices. 

Definition 2.3 (Subsolutions). Let e G Lj oc (R n x (0,T)) with e > 0. A 
subsolution to the incompressible Euler equations with given kinetic energy 
density e is a triple 

(v, u, q) : R n x (0, T) -> R n x 5 nxri x R 

with the following properties: 

• v € Lf oc , u € Lj oc , q is a distribution; 

dtv + div u + V<7 ^ sense of distributions; (6) 
dwv = 0, J > \ j 

v®v — u<^el a.e. . (7) 

Observe that subsolutions automatically satisfy \\v\ 2 < e a.e. (the in- 
equality follows from taking the trace in flTJ)). If in addition we have the 
equality sign *\v | 2 = e a.e., then the v component of the subsolution is in 
fact a weak solution of the Euler equations. As mentioned above, in pass- 
ing to weak limits (or when considering any other averaging process), the 
high-frequency oscillations in the velocity are responsible for the appear- 
ance of a non-trivial Reynolds stress. Equivalently stated, this phenomenon 
is responsible for the inequality sign in ([7]). 

The key point in our approach to prove Theorem 12.21 is that, starting 
from a subsolution, an appropriate iteration process reintroduces the high- 
frequency oscillations. In the limit of this process one obtains weak solutions. 
However, since the oscillations are reintroced in a very non-unique way, in 
fact this generates several solutions from the same subsolution. A brief 
outline of the relevant iteration scheme will be given in Section HI In the 
next theorem we give a precise formulation of the previous discussion. 

Theorem 2.4 (Subsolution criterion). Let e G C(M n x (0, T)) and (v,u,q) 
a smooth, strict subsolution, i.e. 

(v,u,q) e C°°(R n x (0,T)) satisfies © (8) 

and 

v®v-u<^e on R n x(0,T). (9) 
Then there exist infinitely many weak solutions v £ L^ c (R n x (0, T)) of the 
Euler equations such that 
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for a.e. (x,t). Infinitely many among these belong to C((0, T), L 2 ). If in 
addition 

v(-,t) v (-) in L 2 oc (R n ) ast^O, (10) 
then all the v 's so constructed solve ([T]) . 

This theorem is indeed Proposition 2 of [3lj and from it we derive essen- 
tially all the results concerning the incompressible Euler equations which 
are surveyed in this paper. However, for a couple of cases we will need the 
following slightly stronger statement, which can be proved with the same 
methods (see also the Appendix of [75] for a more general statement): 

Proposition 2.5. The same conclusions of Theorem \2.J\ hold if ([8j) and 
([9]) are replaced by the following weaker assumptions 

(1) (v, u, q) is a continuous subsolution; 

(2) for all (x, t) € R n x (0, T] 

-\v\ 2 <e =>■ v ®v — u < 

(3) The domain {(x,t) : ^\v\ 2 < e} C R n x (0, T) has nonempty interior 
and the boundary of each time-slice has Lebesgue measure. 

Condition (2) amounts to the requirement that in the open subset of 
R n x (0, T) where (v,u,q) is not a solution, it should actually be a strict 
subsolution. 

Sketch proof of Proposition \2.5l The argument is exactly the same as that 
of Proposition 2 of [31] . A close inspection of that argument shows that it 
uses only the continuity of the subsolution. It remains therefore to show 
that the open set £1 x (0, T) in Proposition 2 of [31] can be substituted with 
the more general U := {v®v — u < \e\. In the proof of [31] the assumption 
U = O x (0, T) plays a role only in the construction of the grid of Subsection 
4.5 of [31] . In order to handle the case considered here, it suffices to use the 
same cubical partition of W 1 x R and to perform the perturbations only on 
those closed cubes which are contained in U. □ 

2.3. Global existence of weak solutions. One way to utilize Theorem 
12.41 is to start with an initial datum vq € L^ oc (R n ), and construct a smooth 
triple (v,u,q) € C°°(M n x (0, T)) solving © with initial datum vq in the 
sense of (|10p . We further define 

e(v,u) := | \ max {v ®v-u), 

where A max denotes the largest eigenvalue. It is then obvious that (y,u,q) 
is a smooth, strict subsolution for 

e(x,t) := e(y(x,t),u(x,t)^ + min(t, j) 

Of course there are many ways of constructing such a subsolution, since 
the system ([6]) is underdetermined (there is no evolution equation for u\). 
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This observation is closely related to the well known closure problem in 
turbulence. 

By constructing a subsolution with bounded energy, E. Wiedemann in 
|80| recently obtained the following: 

Corollary 2.6 (Global existence for weak solutions). Let Vq € L 2 (T n ) be a 
solenoidal vectorfield. Then there exist infinitely many global weak solutions 
([T]) with bounded energy, i.e. such that 

E(t) = \ [ \v(x,t)\ 2 dx 

Z Jjn 

is bounded. Moreover E(t) —¥ as t — > oo. 

It should be noted, however, that for these solutions the energy E{t) does 
not converge to ^H^oll^ as t — > 0: the energy will, in general, have an 
instantaneous jump. We will return to this issue in the next section. 

2.4. The coarse-grained flow and measure-valued solutions. Follow- 
ing the idea that a subsolution (v, u, q) represents the averaged (or coarse- 
grained) velocity, stress tensor and pressure respectively, it is natural to 
expect that, for any given subsolution, there is a sequence of weak solutions 
Vk with vy. — 1 v in L 2 oc as k —> oo. Indeed, this is a consequence of the proof 
of Theorem 12.41 one obtains a set of weak solutions which is dense in the 
space of subsolutions in the weak I? topology. For details see Section 01 
This can be made even more precise with the concept of Young measures. 

Let us recall that, given a sequence Vk € L 2 oc (M. n x [0,T)), there exists a 
subsequence (not relabeled) and a parametrized probability measure u x t on 
M n x [0, T), called the associated Young measure, such that 

f( v k)~^ / f dv x j in Lfo c for all bounded continuous /. (11) 

Thus, Young measures record the one-point statistics of oscillations in weakly 
convergent sequences. 

In order to capture both high-frequency oscillations as well as possible 
concentrations for Euler flows, DiPerna and Majda introduced the concept 
of measure-valued solutions. With this generalization, one can make sense 
of the limit (jlip also for test functions / which are not necessarily bounded, 
in particular for the stress tensor fij(v) = v l v J . Here we follow Alibert and 
Bouchitte pQ. For such test-functions the limit (jlip takes the form 

/(w fc ) dxdt A ( / / dv x A dxdt + f f f°° dv™\ \{dxdt) (12) 
KJw 1 ' J yJs™- 1 ' J 

where the convergence is in the sense of measures and := lim^oo 

is the L 2 recession function of /. The triple (u, A, v°°) is made of 

• the oscillation measure u x t which is a parametrized probability mea- 
sure on M n ; 
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• the concentration measure A, which is a nonnegative Radon measure 
on R n x (0,T); 

• the concertration-angle measure u^ t which is a parametrized proba- 
bility measure on S 12-1 . 

Note that for bounded / the formula in (|12p reduces to that in (jlip . because 
f°° = 0. 

Definition 2.7 (Measure- valued solutions). A measure-valued solution of 
the incompressible Euler equations is a triple {v x ,t. A, u^ t ) such that the fol- 
lowing two identities hold for all tp G C^°((0,T) x M n ;IR n ) with div(p = 
and for all tp G C c °°((0,T) x R n ): 

/ / d t <p ■ (£, v x j) + Vip : (£ ® £, z/ Xjt ) dxdt 
Jo Jr* ^ (13) 

+ I [ V^> : (0 <g> 0, i£° ) dX = , 

[ (Z,v x , t )-Vipdxdt = 0. (14) 
Jr™ 

y^i) is a measure-valued solution of ([1]) 

if (fl3|) can 6e replaced by 

dtp ■ (£, ^x,i) + V<y9 : (£ <8> £, v x,t) dxdt 

(15) 

Vv9 : (6 ® 0, i/~ f ) dA + / <p(x, 0) • t; (x) dx = 

JR™ 

/or ip G C c °°([0, T) x R n ; R n ) wi/i div v? = 0. 
In these formulas 

v:= {S,u x>t ) = f Sdu x ,t(£) (16) 



(17) 



is the average (coarse-grained) velocity field, and similarly 



u 1 



[ + / 0®6v£ t (6)\(dxdt). 

Jm n J s n - 1 



is the average stress tensor (recall from subsection 12.21 that by u we denote 
the traceless part of the stress tensor). Note that in general the latter is 
a measure rather than a locally integrable function, because of possible L 2 
concentrations. In [33] DiPerna and Majda showed that any sequence of 
Leray solutions of the Navier-Stokes equations converges in the vanishing 
viscosity limit to a measure-valued solution. 

Measure- valued solutions give rise to subsolutions as in Definition 12.31 
Conversely, given a subsolution (U, u, q) , it is trivial to extend it to a measure- 
valued solution by finding for a.e. (x, t) probability measures v Xj t such that 
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()16p and (|17p are satisfied. In this regard it is important to note that in the 
definition of measure-valued solutions there are no microscopic constraints, 
that is, contraints on the distributions of the probability measures. This is 
very different from other contexts where Young measures have been used, 
such as conservation laws in one space dimension |79[ [32] : in these and simi- 
lar references the Young measures satisfy additional microscopic constraints 
in the form of commutativity relations (for instance as a consequence of the 
div-curl lemma applied to the generating sequence). 

Measure-valued solutions are a very weak notion, with a huge scope for 
unnatural non-uniqueness. However, the stronger notion of weak solution in 
Definition 12.11 actually exhibits this very same non- uniqueness, as witnessed 
by the following theorem, proved in |77j . 

Theorem 2.8 (Weak solutions as 1-point statistics). Given a measure- 
valued solution (yxjU^j^xt) °f the incompressible Euler equations, there ex- 
ists a sequence of weak solutions Vk with bounded energy such that (112p 
holds. 



3. Energy 

In the previous section we have seen that weak solutions of the Euler 
equations are in general highly non-unique, at least if we interpret weak 
solutions in the usual distributional sense of Definition 12.11 In particular the 
kinetic energy density \\v\ 2 can be prescribed as an independent quantity. It 
is therefore quite remarkable that, despite this high flexibility, the additional 
requirement that the energy E(t) = \ L„ \v \ 2 dx be non-increasing already 
suffices to single out the unique classical solution when it exists. 

Theorem 3.1 (Weak-strong uniqueness). Let v G L°°([0, T),L 2 (M n )) be a 
weak solution of (pQ) with the additional property that Vv + Vv T € L°° . 
Assume that {y,X,v°°) is a measure-valued solution of (pQ) satisfying 

[ [ \i\ 2 dv x>t {i)dx+ f d\ t {x) < f \v \ 2 (x)dx for a.e. t. (18) 

jR n JR n JR n JR n 

Then (v,X,u°°) coincides with v as long as the latter exists, i.e. 

Vx t = S v(Xtt) for a. a. (x,t) € R n x (0,T) and A = on R n x (0,T). 

This theorem recently appeared in [15], building upon ideas of [HI [16], 
where the authors dealt with the energy of measure-valued solutions to the 
Vlasov-Poisson system. More precisely, the proof of [15] yields the following 
information: if v x ^t satisfies (|18p . then 

v(x,t) := £,dv x>t {ii) (= (£, v x ,t)) 

JR" 

is a dissipative solution of (pQ) in the sense of Lions (see [51]). In fact, Lions 
introduced the latter notion to gain back the weak-strong uniqueness while 
retaining the weak compactness properties of the DiPerna-Majda solutions. 
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Theorem 13.11 shows that this can be achieved in the framework of DiPerna 
and Majda by simply adding the natural energy constraint (118p . 



3.1. Admissible weak solutions. It is easy to see that C 1 solutions of the 
incompressible Euler equations satisfy the following identity, which expresses 
the conservation of the kinetic energy in a local form 



M 2 , ( ( M 2 



d t ^j-+dw I I L Y+P)v) =0. (19) 

Integrating (1191) in space we formally get the conservation of the total kinetic 
energy 

^ [ — {x,t)dx = 0. (20) 
dt Jk™ 2 

These identities suggest that the notion of weak solution to (pQ) can be com- 
plemented with several admissibility criteria, which we list here: 

(a) 

\v\ 2 (x,t)dx< / \vo\ 2 (x)dx for a.e. t. 



(b) 

/M 2 (x,^</MV,»)l<fa fora.e. t >, 
(c) If in addition v S Lf oc , then 

M 2 ( ( M 2 



d t ^-+div I (2" ~*~ P J v ) -° 

in the sense of distributions (note that, since — Ap = divdivu(g)t>, the 
product \v\ 2 p is well-defined by the Calderon-Zygmund inequality). 

Condition (c) has been proposed for the first time by Duchon and Robert 
in [33] and it resembles the admissibility criteria which are popular in the 
literature on hyperbolic conservation laws. 

Next, denote by L^(R n ) the space L 2 (]R n ) endowed with the weak topol- 
ogy. We recall that any weak solution of ([I]) can be modified on a set of 
measure zero so to get v G C([0, T), L^R* 1 ) (this is a common feature of 
evolution equations in conservation form; see for instance Theorem 4.1.1 of 
|29| ) . Consequently v has a well-defined trace at every time and the require- 
ments (a) and (b) can therefore be strengthened in the following sense: 

(a') 



f\ V \ 2 (x,t)dx<f\v \\x)\dx 



for every t. 



(b') 



\v\ 2 (x,t)dx< J \v\ 2 (x, s)\ dx for every t > s. 
However, none of these criteria restore the uniqueness of weak solutions. 
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Theorem 3.2 (Non-uniqueness of admissible weak solutions). There exist 
initial data vo € L°° n L 2 for which there are infinitely many bounded solu- 
tions of dl]) which are strongly L 2 -continuous (i.e. v € C([0, oo), L 2 (W 1 ))) 
and satisfy (a'), (b') and (c). 

The conditions (a'), (b') and (c) hold with the equality sign for infinitely 
many of these solutions, whereas for infinitely many other they hold as strict 
inequalities. 

This theorem is from [31]. The second statement generalizes the intricate 
construction of Shnirelman in [69], which produced the first example of a 
weak solution in 3D of ([I]) with strict inequalities in (a) and (b). 

3.2. Wild initial data. The initial data vn as in Theorem l3.2l are obviously 
not regular, since for regular initial data the local existence theorems and 
the weak-strong uniqueness (Theorem 13. ip ensure local uniqueness under 
the very mild condition (a). One might therefore ask how large is the set 
of these "wild" initial data. A consequence of our methods is the following 
density theorem (cp. Theorem 2 in [77]). 

Theorem 3.3 (Density of wild initial data). The set of initial data vo for 
which the conclusions of Theorem 1 3. ,21 holds is dense in the space of 1? 
solenoidal vectorfields. 

Another surprising corollary is that the usual shear flow is a "wild initial 
data". More precisely, consider the following solenoidal vector field in R 2 

, f (1,0) if x 2 >0 . . 

V °W = { (-1,0) ifx 2 <0. (21) 

or the following solenoidal vector field in T 2 = S 1 xS 1 : 

„ M-J (1.0) if02€(-7T,O) (on) 
M*)- j ( _ 1)Q) if 2e(O)7r ). ( 22 ) 

Theorem 3.4 (The vortex-sheet is wild). For vq as in (|2ip . there are 
infinitely many weak solutions of ([I]) onM. 2 x [0, oo) which satisfy (c). For 
vq as in (I22p there are infinitely many weak solutions of (pQ) on T 2 x [0, oo) 
which satisfy (c), (a') and (b'). 

Theorem 13.41 is proved in [76] using Proposition 12.51 and hence the proof 
essentially amounts to showing the existence of a suitable subsolution. The 
construction of such subsolution follows an idea introduced in [75] for the 
incompressible porous media equation, see Theorem 14.81 below. 

Since the various additional requirements discussed above do not ensure 
uniqueness of the solution even for this very natural initial condition, The- 
orem 13.41 raises the following natural question: is there a way to single out 
a unique, physically relevant solution? In two space dimensions one could 
further impose that the vorticity is a measure, leading to the well-known 
problem of uniqueness for the vortex sheet (we note in passing that our 
methods do not seem to apply to the vorticity formulation of 2-d Euler: 
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compare with the discussion in Section I4.2f) . Two other popular criteria 
considered in the literature for hyperbolic conservation laws are 

(1) the vanishing viscosity limit, 

(2) the maximally dissipative solution. 

For scalar conservation laws they both single out the unique entropy so- 
lution, see [28} 15]. For the 2-dimensional incompressible Euler system the 
situation is surely more complicated. As it happens with Theorem [321 some 
of the solutions constructed in [76] preserve the energy, whereas some other 
are dissipative. On the other hand, it is easy to see that the vanishing 
viscosity limit for initial data as in (I2ip or (|22h is the stationary solution, 
which is obviously conservative. Therefore, even if the two criteria singled 
out unique weak solutions, they would be two different ones. 

In [76] it is also shown that, for the 2-dimensional torus T 2 , the maximal 
dissipation rate reachable with the proof of Theorem 13.41 is 

max f = -1/6. 

It is, however, not clear whether there is a solution with this precise dissi- 
pation rate and, if it exists, whether it is unique. 

4. Applications to other equations 

The ideas introduced in the previous sections apply to many other nondis- 
sipative systems of evolutionary partial differential equations. We start with 
some general considerations, and refer the reader for more details on this 
general framework to [461 156] . 

Following Tartar |79j . we consider general systems in a domain T> C M. d 
of the form 

d 

AdiZ = in V (23) 

i=i 

z(y) G K a.e. y G V (24) 

where 

• z : T> C M. d — > M> N is the unknown state variable, 

• Ai are constant m x N matrices 

• and K C M. N is a closed set. 

Plane waves are solutions of (j23|) of the form 

z(x) = ah(x ■(,), (25) 

where h : R — > R. The wave cone A is given by the states a € such that 
for any choice of the profile h the function ([25]) solves ([23]) . that is, 

A := ja G R N : 3£ G R d \ {0} with J^&^a = j . (26) 

The oscillatory behavior of solutions to the nonlinear problem is then deter- 
mined by the compatibility of the set K with the cone A. This compatibility 
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is expressed in terms of a suitable concept of A-convex hull K A (for the pre- 
cise definition, see Section [473]) . Modulo technical details, the subsolutions 
from Definition 12.31 are solutions z of the linear relations (I23p which satisfy 
the relaxed condition z £ K A . 

The idea of convex integration is to reintroduce oscillations by adding 
suitable localized versions of (|25p to the subsolutions and to recover a so- 
lution of (|23p - (|24p iterating this process. The upshot is that in a Baire 
generic sense, most solutions of the "relaxed system" are actually solutions 
of the original system. 

There are various different forms of implementing convex integration in 
this general framework, see for instance |57l [271 17^1 145] , Common to all 
approaches is that one is working in a space of subsolutions in which highly 
oscillatory perturbations are possible. An elegant way of formalizing this 
was introduced by B. Kirchheim in [J5] Section 3.3. We recall the main 
steps. 

The space of subsolutions arises from a nontrivial open set U C R satis- 
fying the following pertubation property (cp. for instance with Proposition 
2.2 in [30] and the proof of Theorem 4.1 in the same paper). 

Perturbation Property (P): There is a continuous function e : M + — > R + 
with e(0) = with the following property For every z £lA there is a sequence 
of solutions Zj G C^{Bi) of ([23]) such that 

• Z j A in L°°(R d ); 

• z + Zj {y) £U Vy € M. d ; 

• f\zj(y)\ 2 dy >e(dist(z,K)). 

Next, let 

X = {ze C™(V) : ([23]) holds and z(y) € U for all y £ V} 

and let X be the closure of Xq in L°°(2?) with respect to the weak* topology. 
Assuming that K is bounded, the set X is bounded in L°° and the weak* 
topology is therefore metrizable on X. 

An easy covering argument, together with property (P), results in the 
following lemma: 

Lemma 4.1. There is a continuous function e : M + — > R + with e(0) = 
such that, for every z € Xq there is a sequence zj € Xq with 

jT \zj - z\ 2 dy>e (J dist(z(y), K) dy^j 

Since the map z J v \z\ 2 dy is a Baire-1 function on X, an easy applica- 
tion of the Baire category theorem gives that the subset of z € X satisfying 
(|24p is Baire-generic in X. 

The argument just sketched yields weak solutions to (I23p - (l24p (assuming 
that (P) holds for some U) which are zero on the boundary dT> in the trace 
sense with respect to the operator in (l23j) . In particular these weak solutions 
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are extendable by zero to M. d D D. In applications to evolution equations 
T> is a space-time domain, say T> = W 1 x (0,T), and thus this argument 
yields weak solutions with compact time support, as in [66l [26J [71] . For the 
construction of weak solutions with arbitrary initial data, in particular for 
the construction of admissible weak solutions, refinements of this argument 
are necessary. A more detailed exposition for such cases is presented in the 
Appendix of [75] . 

In the following we survey some examples. 

4.1. Compressible Euler. As a byproduct of our analysis of the incom- 
pressible Euler system, non-uniqueness theorems for admissible solutions of 
the so-called p-system were proved in |31j . The system of isentropic gas 
dynamics in Eulerian coordinates is the oldest hyperbolic system of con- 
servation laws and consists of n + 1 equations in n space dimensions. The 
unknowns are the density p and the velocity v of the gas. The equation are 

d t p + div x (pv) = 
d t (pv) + div x (pv ® v) + V[p(p)] = 

p(0,-) = P { 
v(0,-) = v° 

(cf. (3.3.17) in [29] and Section 1.1 of |67| p7). The pressure p is a function 
of p, which is determined from the constitutive thermodynamic relations of 
the gas in question and satisfies the assumption p' > 0. A typical example 
is p(p) = hp" 1 , with constants k > and 7 > 1, which gives the constitutive 
relation for a polytropic gas (cf. (3.3.19) and (3.3.20) of [29]). 

Weak solutions of f)2T[) are bounded functions in W n , which solve it in 
the sense of distributions. Admissible solutions have to satisfy an additional 
inequality, coming from the conservation law for the energy of the system. 

Definition 4.2. A weak solution of (i27|) is a pair (p, v) € L°°(IR n ) such that 
the following identities hold for every test function tp, ip G C£°(IR n x [0, oo[): 

pd t i> + pv ■ V x ifi\ dxdt + / p°(x) tp(x, 0) dx = 0, (28) 
pv ■ dt(f + p(v (g> v,V(f) dxdt + / p°(x)v°(x) ■ (p(x,0) dx =0. 



(27) 



(29) 

Consider the energy e : M + — >• ]R given through the law p(r) = r 2 e'(r). 
A weak solutions of (I27p is admissible if the following inequality holds for 
every nonnegative tp 6 C£°(M. n x M); 



1 1 9 \ / 1 1 2 

p\v\ \ I , . p\v\ 



P<P) + M + ^pe{ P ) + +p(p)j v ■ 

+ / (A(P°) + ^^W,0) > 0. (30) 



The following nonuniqueness result was proved in [31 
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Theorem 4.3 (Non-uniqueness for the p-system). Let n > 2. Then, for any 
given function p, there exist bounded initial data (p } v ) with p° > c > for 
which there are infinitely many bounded admissible solutions (p, v) of (I27p 
with p > c > 0. 

A variant of this theorem has been recently shown by E. Chiodaroli using 
the same techniques (see [20J). Chiodaroli's Theorem highlights that the 
main role in this loss of uniqueness is played by the velocity field. 

Theorem 4.4 (Non-uniqueness with arbitrary density). For every periodic 
p° 6 C 1 with p > c > there exists an initial velocity v° G L°° and a time 
T > such that there are infinitely many bounded admissible solutions (p, v) 
of (]27p on W l x [0, T[, all with density bounded away from 0. 

4.2. Active scalar equations. Active scalar equations are a class of sys- 
tems of evolutionary partial differential equations in n space dimensions. 
The unknowns are the "active" scalar function 6 and the velocity v, which, 
for simplicity, is a divergence-free vector field. The equations are 

8 t o + vve = o 

aiv x v = v ' 

and v and are coupled by an integral operator, namely 

v = T[0\. (32) 

Several systems of partial differential equations in fluid dynamics fall into 
this class. 

We rewrite (f3T1) and ([32]) , in the spirit of Section 12.21 (see also ([6]) and 
0), as the system of linear relations 

dtO + div x q = 

div x v = (33) 
v = T[6] 

coupled with the nonlinear constraint 

q = 6v. (34) 

The initial value problem for the system (|33p -(|34 p amounts to prescribing 
6(x,0) = e (x). 

As described at the beginning of this section, a key point is that the 
linear relations (|33p admit a large set of plane wave solutions. Note that 
these linear relations are not strictly speaking of the form (|23p and in order 
to define a suitable analog of the plane waves in this setting we assume that 
the linear operator T is translation invariant. Let m(£) be its corresponding 
Fourier multiplier. Then we require in addition that 

m(£) is 0-homogeneous (35) 

so that (|33p has the same scaling invariance as (]23p . Furthermore the con- 
straint div^f = implies that 

£ • m(£) = 0. (36) 
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An important remark at this point is that the O-homogeneity of m ex- 
cludes the vorticity formulation of the 2-dimensional incompressible Euler 
equations: indeed convex integration does not seem to apply to this situa- 
tion because the highest order derivatives, namely the vorticity itself, appear 
linearly in the equation. In a geometric context this issue has been raised 
by Gromov in Section 2.4.12 of [40] and to a certain extent analysed in [73j . 
Instead, the 0-homogenity ensures that 9 and v are of the same order. 

In spite of this restriction, several interesting equations fall into this cat- 
egory. Perhaps the best known examples are the surface quasi geostrophic 
and the incompressible porous medium equations, corresponding, respec- 
tively to 

m(0 = ^rV&.Ci) and (37) 
m(0 = |£r 2 (£i6,-£i 2 ). (38) 

In [26] Cordoba, Faraco and Gancedo proved 

Theorem 4.5. Assume m is given by (138|) , Then there exist infinitely many 
weak solutions of ([33]) and ([34]) in L°°(T 2 x [0, +oo[) with 9 = 0. 

This was generalized by Shvydkoy in [71] to all even m(£) satisfying a 
mild additional regularity assumption, namely 

Theorem 4.6. Assume m is even, homogenous and the set 

{?n(£) : m|§n-i is a regular immersion around £} 

spans W 1 . Then there exist infinitely many weak solutions of (133p and (|34j) 
in L°°(T 2 x [0, +oo[) with p = 0. 



4.3. Laminates and hulls. The proofs of Theorems 14.51 and 14.61 in [26\ I71j 

use some refined tools which were developed in the theory of laminates and 
differential inclusions and they present some substantial differences with the 
proofs in [30\ [31] . In order to address these differences we start by recalling 
some of the standard notions in the theory of differential inclusions. These 
notions have been developed in the particular case of gradient vector fields 
(i.e. when (j23f) is given by curl z = 0). Nevertheless the extension to systems 
of the form (I23p - (|24p is straightforward. 

Remark 4.7. In the case of Theorem \4-6\ several modifications of this gen- 
eral strategy are necessary. First of all, as mentioned above, (123j) needs to 
be replaced by (|33p . This in turn means that in general it is not possible to 
obtain a sequence {zj} as in (P) if we insist that each function is compactly 
supported. It is however possible to build up such a sequence if instead of 
requiring that its support be compact, we require that Zj converges uniformly 
to in the complement of -Bi(O) (cp. Lemma 2.1 in [71] ). 
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Returning to the general considerations at the beginning of this chapter, 
we recall that a central point is to find an open set IA satisfying the pertur- 
bation property (P). One possible candidate would be to take the largest 
open set U max satisfying (P). Obviously this set is particularly meaningful 
since it gives the largest possible space X for which the genericity conclusion 
above holds. Moreover, this has the advantage that - at least in many rel- 
evant cases - the set U max coincides with the interior of the lambda-convex 
hull K^, which in turn can be characterized by separation arguments. For 
instance, in Theorem 12.41 condition ([9]) characterizes precisely the interior of 
K A . Moreover, in this case of the interior of is the interior of the convex 
hull K co . As a side remark observe that, by Jensen's inequality, any set U 
with property (P) must necessarily be a subset of the convex hull of K. 

In [26] and [7TJ the authors avoid calculating the full hull and instead 
restrict themselves to exhibiting a nontrivial (but possibly much smaller) 
open set U satisfying (P). However, in exchange they are forced to use much 
more complicated sequences Zj. Indeed, the zj's used in [30] are localizations 
of simple plane waves, whereas the ones used in [26] and |71] arise as an 
infinite nested sequence of repeated plane waves. 

The obvious advantage of the method introduced in [22] and used in |71] 
is that it seems to be fairly robust and general. This is useful in cases where 
an explicit computation of the hull iT A (or even of the convex hull!) is out 
of reach due to the high complexity and high dimensionality. 

On the other hand, the advantage of computing the hull is that it gives 
optimal criteria for wild initial data, as in Section [3.21 Indeed, observe that 
in Theorem 12.41 the initial data vq is given by the subsolution through (|10p . 
In the case of Theorem 14.51 that is, for the incompressible porous medium 
equation, the set K A was computed in [75]. As a consequence one obtains 
the following existence theorem for the Muskat problem (see Theorem 1.1 



Theorem 4.8. Assume m is given by (|38p . There exist infinitely many 
weak solutions of (f33|) and (j34^) in L°°(1R 2 x [0, +oo)) with \6\ = 1 a.e. and 



In fact the solutions constructed in this theorem exhibit a mixing zone (i.e. 
where 6 oscillates wildly between ±1) around the initial interface {x2 = 0}, 
which is expanding linearly in time. 

We recall that the coarse-grained flow for this problem has been analysed 
in detail in|6H 162] , In [61] F. Otto introduced a relaxation approach for 
the incompressible porous media equation, based on the gradient flow for- 
mulation and using ideas from mass transport. It was shown that, under 
certain assumptions, there exists a unique "relaxed" solution 6, representing 
a kind of coarse-grained density. Moreover, Otto proved that, in general, 
the mixing zone (where the coarse-grained density 8 is strictly between ±1) 
grows linearly in time, with the possible exception of a small set of volume 



in [75]). 
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fraction 0{t~ 1 / 2 ) (cp. with [621 Remark 2.1]). Following the general consid- 
erations in Section 12.21 we can interpret 9 as the subsolution corresponding 
to the solutions in Theorem 14.81 f° r which there exists a sequence pk of weak 
solutions such that pk — ~p. 

Although subsolutions are clearly not unique, the coarse-grained density 
of Otto turns out to be extremal in the sense that it corresponds to the 
maximal expansion of the mixing zone. It is interesting to note in this 
connection, that, although weak solutions are clearly not unique, there is 
a way to identify a selection criterion among subsolutions which leads to 
uniqueness. 

5. Gromov's h-principle and the Nash-Kuiper Theorem 

The origin of convex integration lies in the famous Nash-Kuiper theorem. 
In this section we briefly recall some landmark results from the theory of 
isometric embeddings. 

Let M n be a smooth compact manifold of dimension n > 2, equipped 
with a Riemannian metric g. An isometric immersion of (M n ,g) into R m 
is a map u € C 1 (M n ;W m ) such that the induced metric agrees with g. In 
local coordinates this amounts to the system 

diu ■ dju = gij (39) 

consisting of n(n + l)/2 equations in m unknowns. If in addition u is injec- 
tive, it is an isometric embedding. Assume for the moment that g € C°°. 
The two classical theorems concerning the solvability of this system are: 

(A) if m > (n + 2)(n + 3)/2, then any short embedding can be uni- 
formly approximated by isometric embeddings of class C°° (Nash 
[59] , Gromov [40]); 

(B) if m > n + 1, then any short embedding can be uniformly approxi- 
mated by isometric embeddings of class C 1 (Nash |58j . Kuiper [48] ) . 

Recall that a short embedding is an injective map u : M n — > M m such that 
the metric induced on M by u is shorter than g. In coordinates this means 
that 

(d iU ■ d jU ) < ( 9ij ) (40) 

in the sense of quadratic forms. Thus, (A) and (B) are not merely existence 
theorems, they show that there exists a huge (essentially C°-dense) set of 
solutions. This type of abundance of solutions is a central aspect of Gromov's 
/i-principle, for which the isometric embedding problem is a primary example 
(see [401 [36]). 

There is a clear formal analogy between (|39l) - (j40l) and {I]),©, ([7]). First 
of all, note that the Reynold stress measures the defect to being a solution 
of the Euler equations and it is in general a nonnegative symmetric tensor, 
whereas gij — d^u ■ djU measures the defect to being isometric and, for a 
short map, is also a nonnegative symmetric tensor. More precisely (|39|) can 
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be formulated for the deformation gradient A := Du as the coupling of the 
linear constraint 

curl ,4 = 

with the nonlinear relation 

A* A = g. 

In this sense short maps are "subsolutions" to the isometric embedding 
problems in the spirit of Definition 12.31 Along this line of thought, Theorem 
12.41 is then the analogue for the Euler equations of the Nash-Kuiper result 
(B). However note that, strictly speaking, the formal analog of statement 
(B) would be replacing L°° by C° in Theorem 12.41 This type of result is not 
yet available. 

Statement (B) is rather surprising for two reasons. First of all, for n > 3 
and m = n + 1, the system (j39[) is overdetermined. Moreover, for n = 2 
we can compare (B) to the classical rigidity result concerning the Weyl 
problem: if (S 2 ,g) is a compact Riemannian surface with positive Gauss 
curvature and u 6 C 2 is an isometric immersion into M 3 , then u is uniquely 
determined up to a rigid motion ((22j |42], see also |72j Chapter 12 for a 
thorough discussion). Thus it is clear that isometric immersions have a 
completely different qualitative behaviour at low and high regularity (i.e. 
below and above C 2 ). 

A strikingly similar phenomenon holds for the Euler equations: the when 
coupled with the energy constraint \v\ 2 = 2e they are also formally overde- 
termined. Moreover C 1 solutions of the Cauchy problem are unique. We 
will return to these similarities in Section 16.21 

6. Onsager's Conjecture 

The relevance of weak solutions of the Euler equations to turbulence has 
long been surmised. One key question is related to the phenomenon of 
anomalous dissipation. This experimentally observed fact, namely that the 
rate of energy dissipation in the vanishing viscosity limit stays above a cer- 
tain non-zero constant, is expected to arise from a mechanism of transport- 
ing energy from large to small scales (known as an energy cascade) via the 
nonlinear transport term in the Navier-Stokes equations, rather than the 
(dissipative) viscosity term. Motivated by this idea, Onsager conjectured in 
the 1940s that there exist weak solutions of the Euler equations which dissi- 
pate energy. More precisely, the conjecture was that there exist dissipative 
weak solutions in the Holder space C a with a < 1/3, whereas if a > 1/3, 
then the energy is conserved [60l EH] . 

The latter has been rigorously proved with various sharper versions in 
|37[ l24l [34] , in particular we refer the reader to the recent survey [70] . 

On the other hand there is no known construction of a dissipative weak 
solution in C a with any a > 0. As pointed out in Theorem 13.2} one can use 
Theorem l2.4l to produce dissipative solutions which are merely L°°, and prior 
to that the only known construction of a weak solution with dissipation, due 
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to Shnirelman [69] produced an L 2 solution. Despite this, it is instructive to 
take a second look at the construction in Theorem l2.4l in light of expectations 
regarding the energy spectrum and the conjecture of Onsager. 

In Section[4]we presented the so called Baire category method, which is in 
some sense not constructive. However, the same idea of adding oscillatory 
perturbations can be implemented in a constructive way as well. See for 
instance Section 5 of [30]. In a nutshell the idea is to define a sequence of 
subsolutions (v k ,u k ,q k ) as 

(v k+1 ,u k+ i,q k+1 )(x,t) = (v k ,u k ,q k )(x,t) + (V kl U k ,Q k )(x,t,\ k x,\ k t), (41) 

where 

(V k ,U k> Q k )(x,t,£,T) 

is a periodic plane-wave solution of ([6]) in the variables (£, r) with average 0, 
parametrized by (x,t), and X k is a (large) frequency to be chosen. The aim 
is to choose the plane-wave (V k , U k , Q k ) and the frequency X k iteratively in 
such a way that 

• ( v k,u k , q k ) continues to satisfy ([6]) (strictly speaking this requires an 
additional corrector term in the scheme (|4ip ); 

• the inequality ([7]) holds 

• and 

2 

v k ®v k -u k ^- -el. (42) 
n 

Observe that, because of the inequality ([7]), it suffices to show the weak 
convergence in (|42p . 

The role of the frequency can be explained as follows: using (I41D we 
expand as 

v k+ i ® v k+ i - u k+ i = v k ® v k - u k + fv k <g> V k + V k ® v k - U k ^j + V k (g) V k . 

For large A/% the term in brackets converges to zero weakly, whereas the term 
V k <g> V k is non-negative. Therefore, if the sequence {A^l^gN converges to oo 
sufficiently fast, then the sequence of tensors {v k ® v k — u k } is monotone, 
with a uniform bound given by ([7]). The strong convergence follows. 

Thus, convergence of this constructive scheme is improved by choosing 
the frequencies X k higher and higher. In terms of the energy spectrum 
this amounts to eddies (the periodic perturbations) being placed in well 
separated frequencies where the interaction between neighboring frequencies 
is negligible. 

On the other hand clearly any (fractional) derivative or Holder norm of 
v k gets worse by such a choice of X k . The best regularity admitted by the 
scheme corresponds to the slowest rate at which the frequencies X k tend to 
infinity whilst still leading to convergence. 
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6.1. A toy example. 

We will consider here a toy model, where the considerations above can be 
demonstrated and which has been suggested by the anonymous referee to 
[30j . Consider the problem of exhibiting functions u : [0, 1] — > R such that 
\u\ = 1. In the context of (|23p -(|24 p this corresponds to K = { — 1,1} and the 
differential constraint being void. The following scheme aims at producing 
such functions. We assume to start with a given function uq : [0, 1] — > R 
and build a sequence with the following iteration scheme: 

u k+ i(x) = u k {x) + \ [l - u 2 k (x)] s(X k x) , 

where s : R — > R is the 1-periodic extension of l(o,i/2] ~~ 1(1/2,1] an d A& > 1 
is a sequence of frequencies still to be fixed. The following assertions are 
straightforward: 

• If sup[ 0) i] \uk\ < 1, then also supr u |itfc+i| < 1. 

• If supr 0) i] |ito| < 1 and u k — > in L 1 (0, 1), then \uoo\ = 1 a.e. 

Therefore, in order to produce a solution to our toy problem, it suffices to 
choose the sequence {X k } so to ensure the strong convergence of u k . To this 
end observe that 

u k+1 \ 2 dx= / (|ti fc | 2 + \(1 - u\) 2 + u k (l - u 2 k )s(X k x)) dx. 
Jo 

Moreover, as A — > 00, we have s(Xx) — in L 2 (0, 1). Therefore, by choosing 
Afc sufficiently large (depending on u k ), we can ensure that 

/ \u k+ \\ 2 dx > I \u k \ 2 dx + I ^(l — u1) 2 dx. 
Jo Jo Jo 

The strong convergence follows then easily. Here we see that choosing A^ 
to be a rapidly increasing sequence "helps" the strong convergence of the 
scheme. 

However, it is also clear that for any additional regularity of the limit 
one should choose X k to increase as slowly as possible. More precisely, the 
optimal regularity that is reachable via this iteration scheme will depend on 
the connection between the choice of A& with the rate of convergence of the 
scheme. To see this, observe that - roughly speaking - fractional Sobolev 
regularity of will follow from interpolating between the norms 

\\ u k+i — UkWb 1 
\\ u k+l — u k\\BV 



(1 - u\)dx 



Afc / (1 - u 2 k ) dx 



Therefore the following statement is of interest, showing that exponential 
growth of the frequencies leads to exponential convergence of the scheme: 
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Lemma 6.1. Let \ k = 2 k . Then 

\l-u\)dx<{l) k ! (l-ug)ds. 
Jo 

Proof. By the choice of the oscillatory function s we see that for any con- 
tinuous function / E C(— 1, 1) we have 

/ f(u k+1 )dx= [ \[f(u k -\{l-u 2 k ))+\f(u k + \{l-ul))]dx 
Jo Jo 

Now set f(u) := (1 — u 2 ) 1 / 2 . By direct calculation we obtain 

f"(u) = -(l-u 2 )-^ 
and < on (-1, 1). It follows that f"(u)(l - u 2 ) 2 = -f(u) and 

i/(u - «) + \f{u + v)< /(«) + ir(n) V 2 . 
In particular, setting u = ~(1 — -u 2 ) we obtain 

\f(u - |(1 - n 2 )) + §/(« + J(l - u 2 )) < 7 -f(u). 

We conclude 

/ f(u k+1 ) dx < - [ f(u k ) dx 
Jo ° Jo 

and the lemma follows. □ 



6.2. C 1,Q isometric immersions. 

The question of a sharp regularity threshold has been the object of in- 
vestigation for the isometric embedding of surfaces as well (see for instance 
[40] . |81j). As already mentioned, the isometric embeddings of S 2 into M 3 
are rigid in the class C 2 , whereas the h-principle holds for C . Borisov inves- 
tigated embeddings of class C 1,a , and proved the rigidity for a > 2/3 ([8], 
[10] ) and the h-principle for a < 1/13 (although the latter was announced in 
1965, seepT], a partial proof only appeared in 2004 [12]). In [25] we returned 
to this problem, and gave a more modern PDE proof of the h-principle for 
a < 1/7, namely 

Theorem 6.2 (Local existence). Let n E N and go E sym + . There exists 
r > such that the following holds for any smooth bounded open set f2 C M. n 
and any Riemannian metric g E C^(f2) with (3 > and \\g — go\\c° — r - 
There exists a constant 5q > such that, if u E C 2 (Q;W n+1 ) and a satisfy 

i p 



\\ u e ~ 9\\o < an d < a < min 

then there exists a map v E C 1,a (Q; W n+1 ) with 



1 + 2n* 2 



f"e = g and ||« — u||cfi < C ||u"e — g\\^o 
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Corollary 6.3 (Local h-principle) . Let n,go,£l,g,a be as in Theorem \6.2l 
Given any short map u € C 1 (Q;M. n+1 ) and any e > there exists an iso- 
metric immersion v G C 1,a (Q; M n+1 ) with \\u — v\\qo < e. 

The proof of Theorem 16.21 is based on an iteration scheme which follows 
the method of Nash and Kuiper |58l 148] . The iteration consists of stages, 
and each stage consists of several steps. The purpose of a stage is to correct 
the error g — u^e. In order to achieve this correction, the error is decomposed 
into a sum of primitive metrics as 

ii,„ 

g - u i e = ^2 a\v k <g> v k . (43) 
fc=i 

To keep the notation simpler, in what follows we will use the abbreviations 
|| • ||tv and || • \\k, a for the C N and C ,a norms respectively. The natural 
estimates associated with (f43l) are 

ii n n it H 1/2 

l|ofc||o ~ IIS - u e llo 

IKHat+i ~ IMk+2 for N = 0,1,2, 

A step then involves adding one primitive metric. In other words the goal 
of a step is the metric change 

u"e h4 u^e + a 2 v ® v. 

Nash used spiralling perturbations (also known as the "Nash twist") to 
achieve this; in codimension one Kuiper replaced the spirals by corrugations: 
in both cases the new map u is of the form 

u(x) = u(x) + —w(x, Xx ■ v) 

A 

for some appropriate choice of w, where w is periodic in the second variable. 
The parameter A is given by 

X = K- N|2 



where K is a large (but fixed) constant. One checks that the addition of a 
primitive metric is then possible with the following estimates: 

C°-error in the metric ~ \\q — it"e||n — 

" 11 K 

increase of C 1 -norm of u ~ \\g — u"e||o 

increase of C 2 -norm of u ~ 1 1 1 1 2 K 

for any K > 1. Observe that the first two of these estimates is essentially 
the same as in |58U48j . Furthermore, the third estimate is only valid modulo 
a "loss of derivative". The exponent 1/2 in the second estimate is due to 
the quadratic nature of the nonlinearity. 

The low codimension forces the steps to be performed serially. Thus the 
number of steps in a stage equals the number of primitive metrics in the 
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above decomposition which interact. This is given by = n{n + l)/2. To 
deal with the " loss of derivative" problem in |25] we mollify the map u at the 
start of every stage, in a similar manner as is done in a Nash-Moser iteration 
(see |54j . [55j). A careful estimate (see Proposition I6.4j) shows that this can 
be done without any additional error. Therefore, iterating the estimates for 
one step over a single stage (that is, over n* steps) leads to 

C°-error in the metric ~ ||<? — u"e||o — (44) 

K 

increase of C 1 -norm of u ~ \\g — v}e\\^ 2 (45) 

increase of C 2 -norm of u ~ IMI2 K n * (46) 

With these estimates, iterating over the stages leads to exponential conver- 
gence of the metric error, to a controlled growth of the C 1 norm and to an 
exponential growth of the C 2 norm of the map. In particular, interpolation 
between these two norms leads to convergence in C 1,Q for a < Trk — = 

° l+2n* 

1 

l+re(n+l) ■ 

Recalling the analogy with the Euler equations, which amounts to identi- 
flng the velocity v with the deformation gradient Du, the analogous scheme 
for producing C 0,a solutions of Euler would involve perturbations of the 
form 

v(x, t) = v(x, t) + w(x, t, Ax, At) 

(with w periodic and average in the second variable). If the estimates 
analogous to (l4"4"l) - ([4*6]) would carry over with one single step in a stage (i.e. 
n t = 1), one would reach the exponent conjectured by Onsager. This ansatz 
for the critical exponent 1 /3 is independent of considerations on the energy 
spectrum, giving an entirely new point of view on the famous Kolmogorov 
5/3 law. 

6.3. Conservation of energy and rigidity. 

We turn next to the Weyl problem, i.e. the characterization of isometric 
embeddings of 2-dimensional positively curved surfaces in M 3 . Concerning 
the rigidity, it is known from the work of Pogorelov and Sabitov that 

(1) closed C 1 surfaces with positive Gauss curvature and bounded ex- 
trinsic curvature are convex (see |64J); 

(2) closed convex surfaces are rigid in the sense that isometric immer- 
sions are unique up to rigid motion [63] ; 

(3) a convex surface with metric g € C 1 '" with I > 2, < (3 < 1 and 
positive curvature is of class C 1,13 (see [MIES])- 

Thus, extending the rigidity in the Weyl problem to C 1,a isometric immer- 
sions can be reduced to showing that the image of the surface has bounded 
extrinsic curvature in the sense of Pogorelov. This concept is easily ex- 
plained in the following way. Let (M, g) be a smooth 2-dimensional closed 
Riemannian manifold and u : M — s> R 3 an isometric immersion. According 
to Gauss' original definition, the Gauss curvature of u(M) can be defined 
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as the area distortion of the Gauss map N. In a weak sense, this is a well- 
defined object even when the immersion is merely a C 1 map, because in 
this case iV is continuous. The immersed surface u(M) is then said to have 
bounded extrinsic curvature if this area distortion is bounded. 

Using geometric arguments, in a series of papers [H El El El [10] Borisov 
proved that for a > 2/3 the image of surfaces with positive Gauss curvature 
has indeed bounded extrinsic curvature. Consequently, rigidity holds in 
this range and in particular 2/3 is an upper bound on the range of Holder 
exponents that can be reached using convex integration. 

A short proof of Borisov's rigidity result was provided in [25]. Note that 
if u G C 3 one can compute the area distorion of the Gauss map from the 
Riemann-curvature tensor, which in turn depends only on the metric: we 
call this scalar field k. Even if the metric g is smooth, nonetheless this 
identity is in general false if the isometry is not regular enough, as shown 
precisely by the Nash-Kuiper theorem. However, in [25] we showed that 
the C 1, 3 +e regularity is enough for this identity to hold. A key observation 
of the proof in [25] is that this identity can be expressed in the following 
integral formula 

/ f(N(x)) K (x)dA(x) = I f(y)deg(y,V,N)da(y) (47) 
Jv Js 2 

where 

• V is an arbitrary open subset of M; 

• / is any C 1 function on S 2 ; 

• deg (y, V, N) is the Brouwer degree of the map N\y at y. 

The identity (|47p is just a change of variable formula and, surprisingly, 
it can be proved for immersion u G C ,1,2 / 3+e by a regularization argument. 
The corresponding computations in [25] rely on commutator estimates which 
are surprisingly close to the argument of [M] for proving the energy conser- 
vation of C 0,a solutions of Euler when a > 1/3. In the case of isometric 
embeddings there does not seem to be a universally accepted critical expo- 
nent (cp. Problem 27 of [81]), even though 1/2 and 1/3 both seem relevant 
(cp. with our discussion below and with the discussion in [E]). The key 
commutator estimates of |25j are given by the following proposition, which 
can be compared to [24] . 

Proposition 6.4 (Quadratic estimate). Let £1 C K™ be an open set, v G 
C 1 ' a (^l,W n ) and e the Euclidean metric on R m . Consider the corresponding 
pull-back metric (v"e)y = diV ■ djV. Assume v^e G C 2 and let ip G C^°(R n ) a 
standard symmetric convolution kernel. Then, for every compact set K C 0, 

\\{v* W )h-vh\\ c , {K) = Oi^' 1 ). (48) 

In particular, fix a map u and a kernel tp satisfying the assumptions of 
the Proposition with a > 1/2. Then the Christoffel symbols of (v * tp^e 
converge uniformly to those of ire. This corresponds to the results of Borisov 
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in [HJ [7], and hints at the absence of /i-principle for C '2 +£ immersions. One 
might further notice that the regularity C 1 '^ 6 is still enough to guarantee 
a very weak notion of convergence of the corresponding Christoffel symbols. 

7. Outlook and further problems 

We list in this section several open questions which remain to be ad- 
dressed. 

Problem 1. Analogues of Theorem \3.4\ for compressible Euler. 

Different behaviours are expected for different types of Riemann problems. 
The direct analogue of Theorem 13.41 would be a contact discontinuity, and 
it seems likely that some form of the statement survives for compressible 
Euler. Nevertheless, because of the different role of the pressure term, a 
direct transfer from the incompressible to the compressible case (as was 
done for Theorem 14.31 in |31] ) seems impossible. 

Problem 2. Compute the hull for compressible Euler. 

The proof of Theorem 14.31 essentially relies on the observation that The- 
orem 12.41 gives weak solutions of the incompressible Euler system (TjQ) with 
fixed \v | 2 and fixed pressure p. Thus, an appropriate choice automatically 
yields weak solutions of the compressible Euler system (|27p . Within the gen- 
eral framework of ()23|) - (|24[) this route via the incompressible Euler equations 
amounts to restricting to a smaller constitutive set K for compressible Eu- 
ler, where p is fixed to be constant. A further step was obtained recently 
by Chiodaroli in Theorem 14.41 but the general case, which would allow for 
wild oscillations in the density as well as the velocity, remains open. In 
other words, one should define the state space and the constitutive set for 
compressible Euler in full generality, and compute the hull. This would lead 
to the correct definition of subsolution, and thus eventually to the solution 
of the next problem: 

Problem 3. Analogue of Theorem \2. <§1 for compressible Euler. 

As noted in Section 16.31 although there is a surprising number of simi- 
larities between weak solutions of the Euler equations and the Nash-Kuiper 
theorem, in contrast to Euler there are no natural conjectures for critical 
Holder exponents for isometric embeddings. Therefore, a central question 
is: 

Problem 4. Find the optimal Holder- exponent for Theorem 1 6. £1 and Corol- 
lary E3 

It should be pointed out that for incompressible Euler the exponent 1/3 
is conjectured to be critical with respect to the energy conservation. It is 
not known whether being above this exponent implies uniqueness. On the 
other hand the rigidity statement in Section [6.31 is a uniqueness statement. 
A conceivable scenario is then the existence of other regimes between the 
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exponents 1/3 and 2/3 for isometric embeddings, where the rigidity does 
not hold but the /i-principle is also excluded. 

In Section 16.21 we restricted our attention to isometric embeddings in 
codimension 1. Looking closely at the construction it is natural to expect 
that higher codimension for the embedding should result in better regularity. 
Indeed, for the case of sufficiently large codimension A. Kallen showed in 
|43j that the regularity in Theorem 16.21 and Corollary 16.31 can be pushed up 
to any a < 1, provided the metric is sufficiently smooth! 

Problem 5. Analyze the effect of higher codimension on the regularity in 
Theorem 1 6. 2\ and Corollary \6. 31 

In other words, interpolate between Theorem 16.21 (codimension 1) and 
Kallen's result (codimension sufficiently large). 

Regarding Onsager's conjecture the next important question is how to 
modify convex integration so that the weak solutions in Theorem 12.41 are 
Holder continuous. 

Problem 6. Construct continuous and Holder- continuous solutions of in- 
compressible Euler. 

An explicit construction could lead to insights concerning several assump- 
tions usually made in the statistical theory of turbulence (such as the local- 
ness of interactions between eddies, cp. [19], or the very idea of an energy 
cascade). It would also open the way towards proving a true homotopy- 
principle for weak solutions, in the long run perhaps shedding light on the 
global topology of the moduli space of solutions. 

Problem 7. Calculate /estimate the maximal dissipation rate for wild initial 
data. 

By a simple contradiction argument it can be seen that there has to be a 
finite maximal rate. A first example might be the shear flow, cp. with Theo- 
rem [331 For the very similar problem of the incompressible porous medium 
equation, see Theorem 14.81 the maximal expansion rate of the mixing zone 
has been obtained in [62J. 

Apart from a dissipation rate principle, the natural candidate for a selec- 
tion principle for the incompressible Euler equations is the vanishing viscos- 
ity limit. Of course this has several major technical stumbling blocks, but 
already a "recovery sequence" in the vanishing viscosity limit would be of 
interest: 

Problem 8. Recovery sequence for weak solutions of Euler from Navier- 
Stokes. 

It is conceivable, that one would need to consider Navier-Stokes with a 
random perturbation. Indeed, adding noise is not just physically motivated, 
but in some sense necessary: without noise the statement is false (see The- 
orem 13.411 . 
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Problem 9. Probabilistic approach to convex integration. 

Convex integration can be seen as a control problem: at each step of the 
iteration, one has to choose an admissible perturbation, consisting essentially 
of a (plane-)wave direction and a frequency. The following is therefore a 
natural question: are there meaningful optimality conditions? 

Finally, up to now there has been no attempt to link weak solutions and/or 
subsolutions to generalized flows or to Brenier's sharp measure- valued solu- 
tion p3]: 

Problem 10. Relation of subsolutions to generalized flows. 
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